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It is known that, after formation, a Rydberg wave packet undergoes a 
series of collapses and revivals within a time period called the revival 
time, t rev , at the end of which it is close to its original shape. We study 
the behavior of Rydberg wave packets on time scales much greater than 
trev We show that after a few revival cycles the wave packet ceases to 
reform at multiples of the revival time. Instead, a new series of collapses 
and revivals commences, culminating after a time period t ST ^> t TCV 
with the formation of a wave packet that more closely resembles the 
initial packet than does the full revival at time t iev . Furthermore, at 
times that are rational fractions of t sr , the square of the autocorrelation 
function exhibits large peaks with periodicities that can be expressed 
as fractions of the revival time t rev . These periodicities indicate a new 
type of fractional revival occurring for times much greater than t Tev . A 
theoretical explanation of these effects is outlined. 
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An electron wave packet in a Rydberg atom provides a physical system in which 
to explore the interface between classical and quantum mechanics. When a Rydberg 
atom is excited by a short laser pulse, a state is created that has classical behavior for 
a limited time 0. Depending on the excitation scheme, either a radially localized 
wave packet in an eigenstate of angular momentum [^, [|, |5|] or a packet localized 
in the angular coordinates |7| is produced. Either a radial or a circular wave 
packet initially oscillates with the classical keplerian period T c \ = 27m 3 , where n is 
the central value of the principal quantum numbers excited in the packet. However, 
after a finite number of classical orbits quantum-interference phenomena appear and 
the wave packet collapses. 

On a time scale large compared to the classical period T cl , a Rydberg wave packet 
passes through a sequence of fractional and full revivals |], p], @> ^ rev i va l 
occurs at a timeQ t rev = %T c i, when a packet close to the original shape reappears. 
The fractional revivals occur earlier than this, at times that are rational fractions of 
the revival time t rev . They correspond to the formation of macroscopically distinct 
subsidiary wave packets that oscillate with periodicity that is a fraction of the clas- 
sical orbital period T c \ ||. Recent experiments have detected fractional revivals with 
periodicities as small as T rs | T c \ [ |T0|| . 

In this letter, we examine the time evolution and revival structure of both radial 
and circular Rydberg wave packets on time scales much greater than the revival time 
t rev . A new system of full and fractional revivals is uncovered, with structure different 
from that of the usual fractional revivals. 

In Ref. [[I]], long-term revivals at t = n 2 T c \ and t = n 3 T c \ were found. They were 
obtained by expanding the energy in a Taylor series in n through finite order and 
determining the times at which the time-dependent phase in the wave function is an 



integer multiple of 2rr. Ref. [1 !| generalized this approach and found a hierarchy of 
recurrence times for full revivals of the wave packet for t 3> t rev . 

Our approach here is different. We do not look for commensur ability of the terms 



1 Note that t ICV is often denned as half the full revival time, since this is the time at which a 
full wave packet reforms for the first time, though with a phase difference relative to the classical 
motion. We adopt the notation of Ref. [|| for labeling the fractional revivals and use their definition 

Of trev 
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in the phase with integer multiples of 2ir. Rather, we show that at certain times 
greater than t rev the wave function can be expanded in a set of subsidiary waves. 
This method is suitable for studying fractional revivals || as well as full revivals. It 
also reveals the presence of a long-term full revival occurring earlier than those found 



previously and not among the hierarchy of levels noted in Ref. ITT 



A simple way of demonstrating the formation of the new fractional and full revivals 
is to construct the autocorrelation function of the time-evolved wave function with 
the initial packet. At the various fractional- revival times, the absolute square of the 
autocorrelation function exhibits periodicities reflecting the presence of underlying 
subsidiary wave packets. Furthermore, the ionization signal in a pump-probe exper- 
iment should display the same periodicities as the autocorrelation function. Thus, 
the primary features of our results relevant for experiments on Rydberg wave packets 
can be exhibited directly in a graph of the autocorrelation function as a function of 
time. For this reason, we first present plots of the autocorrelation function for times 
greater than £ rev . We use these to indicate the structure of the new revivals. We then 
outline the theoretical explanation of the effects and suggest experiments to detect 
them. 

Provided the eigenenergies are independent of the angular-momentum quantum 
numbers, it can be shown that the autocorrelation function for a given radial packet 
is the same as that for a related circular packet. For example, in the case of hydrogen 
the expression for the modulus squared of the autocorrelation function is 

l^)iH£KI 2 ^T C 1 ) 

n 

for both types of wave packet. Here, c n = (^/(0)\^/ n ), and \I/ n is an energy eigenstate 
with principal quantum number n. 

We consider Rydberg wave packets of hydrogen for the case where the excitation 
spectrum is strongly peaked around a central value n. For simplicity in the present 
work, we take a normalized gaussian distribution of width a for the modulus squared 
of the coefficients, |c n | 2 . However, the results we obtain hold equally well for asym- 
metric excitation distributions, such as those occurring in a description in terms of 



radial squeezed states [O. Our results also can be generalized to the case of nonin- 
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teger values of n, including thereby the effects of quantum defects. This means, for 
example, that our results apply equally well to alkali-metal atoms. Experimentally, 
the value of a is set by the exciting laser pulse. In what follows, the values chosen for 
a correspond to a laser-pulse duration that is effective in displaying the new revival 
structure. 

The time-dependent wave function for the Rydberg wave packet may be written 
as an expansion in terms of eigenstates of hydrogen 

oo 

V(?,t)= c k <p k (f)exp(-iE n t) , (2) 

k=— oo 

where k = n — n, n is assumed large, and <Pkif) is a hydrogenic wave function. Since 
we are taking a distribution strongly centered around n, we may expand E n around 
n: 

oo 

*(r,t) = E c fc <^fc(r)exp 

k=— oo 

In this equation, we have introduced three distinct time scales T c \, t rev , and t sr that 
naturally appear in the first three orders in the expansion. The first two are the 
usual time scales relevant in the description of the conventional revival structure. 
The smallest time scale, T c \, is the classical period for an electron in a keplerian orbit. 
If this term alone appeared in the expansion, the energies would be equally spaced 
and the wave packet would undergo simple harmonic motion without changing shape. 
The second term defines the revival time scale, t rev , and is responsible for the collapse 
and fractional/full revivals of the wave packet ||. 

The third-order term is included because in this paper we are interested in times 
much greater than the revival time. This term defines a new time scale, which we refer 
to as the superrevival time t sr . In terms of the revival time, it is given by t ST = ^t Tev . 
For typical values of n, t sr is between one and two orders of magnitude greater than 
the revival time t IeY . However, this is still several orders of magnitude less than 
the lifetime of the excited Rydberg atom. It is therefore reasonable to examine the 
behavior of the wave packet and determine the effects of the third-order term on the 
revival structure up to times of order £ sr . 

With the gaussian distribution for \c n \ 2 discussed above, we can evaluate the 
absolute square of the autocorrelation function |A(t)| 2 directly from Eq. (|l|) and plot 
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the result. We begin with an example involving a circular wave packet that serves to 
illustrate some key features of the new revival structure. Subsequently, we examine 
a radial wave packet with a lower value of n of a type that could be experimentally 
created and observed. 

Consider the circular wave packet treated in Ref. 0, with n = 320, I = m = n — 1, 
and a = 2.5, which displays a clear sequence of conventional fractional/full revivals 
on the revival timescale t Tev ~ 1.06 //sec. Figure 1 shows the absolute square of 
the autocorrelation function for a Rydberg wave packet with n = 320 for a total 
time period large compared to the revival time. In this case, //sec, and 

we consider times up to just beyond ~£ sr . Figure la shows the behavior of the wave 
packet for approximately the first 9 //sec. The full revival appears at t iev ~ 1.06 
//sec, and fractional revivals are evident at fractions of this value. The largest revival 
appears at t — |t rcv , where the wave function has reformed into a single wave packet, 
but at a time that is out of phase with a classically propagated distribution . 

As can be seen from the figure, there are four or five full revival cycles, with 
peaks gradually decreasing in size and losing the canonical periodicity. However, 
after about 7 //sec a new structure with a noncanonical periodicity is visible in the 
autocorrelation function. Similar structures reappear later, at times near 14, 21, 
28, and 42.5 //sec in Figures lb, lc, Id, and le, respectively. The amplitudes and 
periodicities of these structures change with time. At t ~ 42.5 //sec, a series of peaks 
appears with maximum amplitude greater than that of the full revival at t = t iey . 

The features of this illustrative example can be predicted theoretically. We have 
proved that at times t ~ H ST , where q must be an integer multiple of 3, the wave 
packet can be written as a sum of macroscopically distinct wave packets. The proofs 
are rather technical and are presented elsewhere Hl3 |. Explicitly, we find that near 
times t ~ -t ST we may write 

^(r,t) = YbsMr,t + fT d ) , (4) 

s=0 

where the macroscopically distinct wave packets are defined by 

oo 

M?,t)= J2 c fc y?fc(f)exp [-27TZ (^)] (5) 

k=— oo 



4 



and where their weights b s in the sum (|]) are given by 




) 



(6) 



The quantities a and I are integer constants that depend on q and n: 



2n 




q if q/9 ^ (mod 1) 
q/3 if q/9 = (mod 1) 



(7) 



where N is an integer consisting of the product of all factors of 2n that are also factors 



We have shown that the properties of the b s coefficients are such that the auto- 



a much greater time scale than the usual fractional revivals, which have periodicities 
that are fractions of T c \. Moreover, for special values of q the number of nonzero b s 
coefficients is small and one coefficient dominates. This means that the wave function 
is mostly concentrated in one packet, which explains the peaks in the autocorrelation 
function. In special cases, only one nonzero b s value exists. This first happens when 
t ~ |^sr ; so that q = 6. It corresponds to the formation of a single packet, i.e., a 
full superrevival, with periodicity T pa ^t rev in the autocorrelation function. For the 
fractional times t ~ H sr with q > 6, we find that more than one b s coefficient is 
nonzero, corresponding to the formation of distinct subsidiary wave packets at times 
t 3> trcv We refer to these as fractional superrevivals. 

Since the b s coefficients in (|j) have moduli that are generally not equal, this indi- 
cates that the fractional superrevivals consist of an unequally weighted superposition 
of subsidiary wave packets, which is a feature that is not true for the usual fractional 
revivals. Moreover, after the formation of a fractional superrevival consisting of sev- 
eral subsidiary wave packets, it often happens that the subsidiary packets quickly 
evolve into a configuration where one of them is much greater than the others. The 
dominant wave packet in this case can again resemble the initial wave packet more 
closely than the full revival does at time t rev . In this way, large peaks can form in the 
autocorrelation function near the fractional superrevival times. 

2 The actual periodicity is T = -t Tev + ^T c \, where u and v are integers that depend on n and q. 
Since T c \ <C t Tev , we write T « |t re v 
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correlation function is periodic, with a periodF] T pa -t rev . These periodicities are on 
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We demonstrate these results for the illustrative example with n = 320. In this 
case, the theory predicts the first full superrevival appearing at t ~ 42.5 //sec, with 
an autocorrelation-function periodicity T ~ 530 nsec. This agrees with the structure 
seen in Figure le. For q = 36, 18, 12, and 9, fractional superrevivals are predicted 
to appear, since for these values of q more than one of the b s coefficients is nonzero. 
The corresponding times at which structures should be seen are t ~ 7.08, 14.2, 21.2, 
and 28.3 //sec, and they should have autocorrelation-function periodicities T ps -^t Tev , 
l^rev, jt Tev , and ht Tev , respectively. These predictions agree with the structures seen 
in Figs, la - d. 

To obtain further insight into the formation of a superrevival, we can examine 
the cross-sectional view of a circular wave packet at various times. Figure 2 shows 
a cross-sectional slice of the circular wave packet with n = 320. We view the wave 
packets in the plane of the orbit and plot a circular slice *&(</>) at a fixed radius given 
by the expectation value (r) = |n(2n + 1). Figure 2a shows the initial wave packet 
at time t = 0. Figure 2b shows the wave packet at the canonical revival time t ps t Tev 
when the autocorrelation function is a local maximum. As can be observed, the wave 
packet has indeed reformed and resembles the original wave, but it is asymmetric 
and several smaller subsidiary wave packets also appear. Figure 2c shows the full 
superrevival at t ~ |t sr , at which point only one b s coefficient in Eq. (|4]) is nonzero. 
In comparing Figures 2b and 2c, the superrevival wave packet is seen to resemble more 
closely the initial wave packet of Figure 2a than does the canonical revival packet. 

We chose n = 320 as an illustrative example both for purposes of comparison 
with previous results for t < t rev and because the full and fractional superrevivals are 
more evident for larger n. However, similar effects can be observed for Rydberg wave 
packets with relatively small values of n. 

Figure 3 shows the square of the autocorrelation function for hydrogen with n = 
48 and o = 1.5. As remarked above, the autocorrelation function for hydrogen is 
the same for both radial and circular wave packets. Thus, Figure 3 represents the 
autocorrelation function of a radial Rydberg wave packet that has been excited by 
a short laser pulse centered on the value n = 48. In this case, the full revival is 
at t ps t iev ~ 0.538 nsec. For times greater than this, one can observe a fractional 
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superrevival at t ~ j^t ST ~ 1.61 nsec with autocorrelation periodicity T « |i rev and 
a full superrevival at t ~ |i sr ~ 3.23 nsec with autocorrelation periodicity T « ~t rev . 
The size of the peak in the autocorrelation function shows that the superrevival 
resembles the initial wave packet more closely than does the revival wave packet at 

v ~ trey 

We have proved that similar results hold when quantum defects are present in the 
energies. Therefore, we expect superrevivals to occur for wave packets in alkali-metal 



atoms as well as in hydrogen []14]]. In alkali- metal atoms, the associated characteristic 
times depend on n* = n — 5(1) instead of n, where 5(1) is the quantum defect. 

Figure 3 indicates it is likely an experiment can be performed to detect the full and 
fractional superrevivals discussed in this paper. One possibility is to use the pump- 
probe time-delayed photoionization method of detection for radial Rydberg wave 
packets excited in alkali-metal atoms with n « 45 - 50. Single-photon absorption 
would produce a packet with localization purely in the radial coordinate and with 
the angular structure of a p state. The procedure is experimentally feasible, provided 
a delay line of 3 - 4 nsec is installed in the apparatus. In such an experiment, the 
fractional superrevival at y^tsr, with periodicity T « jt Tev , could be observed with a 
delay line of approximately 1.5 nsec. For even smaller values of n, the required delay 
times can be reduced below 1 nsec. With n ~ 36, for example, the full/fractional 
superrevivals could be detected with delay lines used currently in experiments. 

We enjoyed conversations with Charlie Conover and Duncan Tate. R.B. thanks 
Colby College for a Science Division Grant and the Theory Group at Indiana Univer- 
sity for its hospitality. 
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Fig. 1: The absolute square of the autocorrelation function for a Rydberg wave packet 
with n = 320 and a = 2.5 is plotted as a function of time in microseconds. The 
value a = 2.5 for a wave packet with n = 320 corresponds to excitation with 
a 160 picosecond laser pulse, (a) < t < 9 /isec, (b) 9 < t < 18 /isec, (c) 
18 < t < 27 usee, (d) 27 < t < 36 //sec, (e) 36 < t < 45 ^sec. 

Fig. 2: Unnormalized circular wave packets with n = 320 and a = 2.5. Cross- 
sectional slices of the wave packet in the plane of the orbit and for r = (r) = 
|n(2n + 1) are plotted as a function of the azimuthal angle <p in radians, (a) 
t = 0, (b) t ~ t TCV , (c) t ~ |t sr . 

Fig. 3: The absolute square of the autocorrelation function for a Rydberg wave packet 
with n = 48 and a = 1.5 is plotted as a function of time in nanoseconds. The 
value a = 1.5 for a wave packet with n = 48 corresponds to excitation with a 
900 femtosecond laser pulse. 
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